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We study the algebraic decay of the survival probability in open hierarchical graphs. We present 
a model of a persistent random walk on a hierarchical graph and study the spectral properties of 
the Frobenius-Perron operator. Using a perturbative scheme, we derive the exponent of the classical 
algebraic decay in terms of two parameters of the model. One parameter defines the geometrical 
relation between the length scales on the graph, and the other relates to the probabilities for the 
random walker to go from one level of the hierarchy to another. The scattering resonances of the 
corresponding hierarchical quantum graphs are also studied. The width distribution shows the 
scaling behavior P(T) ~ 1/T. 
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I. INTRODUCTION 

Typical Hamiltonian systems are non-integrable and have a mixed phase space, where regions of regular and chaotic 
motions coexist. The chaotic dynamics of mixed systems is clearly different from the fully chaotic case. This is manifest 
in the behavior of the survival probability in open systems. 

Assume we have an infinite hierarchy of Kolmogorov- Arnold- Moser (KAM) small islands interspersed in a connected 
chaotic region, and suppose we draw a boundary at a given level of this hierarchy, such that the particles leaving 
this boundary are lost. Consider a large initial number No of randomly chosen (with respect to a given probability 
distribution) initial conditions in the chaotic region and let them evolve by the dynamics up to some time t. The 
survival probability P(t) is the ratio N(t) /No in the limit of large No, where N(t) is the number of particles remaining 
within the boundary at time t. In the typical case this probability is believed to decay algebraically, 

t-s 
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It has been argued that the algebraic decay is due to the hierarchical structure of phase space jjj | , 
despite significant efforts, the mathematical understanding of the behavior described by Eq. (1 
Much of our current knowledge of this problem is based on the self-similar Markov chain model 
an expression for the exponent 5 in terms of the parameters of the model. Yet a precise and simple understanding of 
the mechanism based on dynamical properties is lacking. 

In fully chaotic open systems, the survival probability decays exponentially, 

P(t) ~ e-v* . (2) 

This case is well understood. The evolution operator of the probability densities, the Frobenius-Perron operator, 
admits a spectral decomposition in terms of Pollicott-Ruelle resonances 1 || which characterize the relaxation proper- 
ties. In particular, for open systems, the leading resonance is identified as the escape rate 7 in Eq. (^J), and describes 
the slowest relaxation mode of the probability distributions. We point out that in closed systems an equilibrium 
state exists (the leading resonance is equal to zero), and one can study the relaxation to this equilibrium state by 
considering the next leading resonance. In contrast, for open systems the final state does not exist due to the escape, 
the rate of which is characterized by the leading resonance of the Frobenius-Perron operator. We refer to |lCj ] for 
more details concerning the connection between open and closed systems. 

The escape rate can also be interpreted as a macroscopic quantity resulting e. g. from a diffusion process described 
by a Fokker-Planck equation for the macroscopic density of particles. This connection between microscopic dynamics 
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and macroscopic processes, known as the escape rate formalism g |(J 0, Q ||, yields expressions of the transport 
coefficients, e. g. the diffusion coefficient, in terms of the dynamical quantities. The existence of this connection relies 
heavily on the hyperbolic properties of the system, i. e. (i) (almost) every point in phase space is assumed to be of 
saddle type, and (ii), for the open boundaries, the repeller is fractal. 

The absence of an exponential decay rate of the survival probability for a typical system Eq. ([!]) is associated to 
anomalous transport, e. g. in a diffusive process the mean square displacement grows with a power of t not equal to 1. 
With this respect, the connections between macroscopic phenomena and microscopic dynamics are far less understood 
for typical systems than they are in the fully chaotic case. 

Attempts to describe the relaxation properties of systems with mixed phase space in terms of spectral properties of 
the Frobenius-Perron operator have introduced regularization procedures which amount to truncating the Frobenius- 
Perron operator in a finite matrix representation |14[ ]. An alternative approach |L5j considers the presence of a 
vanishing noise and yields finite values of the leading relaxation rates. In both these approaches, it is worthwhile 
stressing that the relaxation rates are the analogues of the leading Pollicott-Rucllc resonances mentioned above in 
reference to relaxation in fully chaotic systems. 

Our purpose in this paper is to understand what properties of the Pollicott-Ruelle spectrum characterize the 
algebraic as opposed to exponential decay of the survival probability. We will do so by considering a model whose 
finite approximations are fully chaotic, but which displays algebraic decay of the survival probability as a limiting 
property. 

For the purpose of this endeavor, we propose to investigate the decay properties of an open onc-dimensional 
hierarchical graph, whose survival probability turns out to decay algebraically. A graph is a collection of bonds on 
which a classical particle has a uniform one-dimensional motion. The bonds are interconnected by vertices where 
neighboring bonds meet. At the vertices, the particle undergoes a conservative collision process with the result that 
its velocity may change direction. In practice this collision process is determined by a random process where outputs 
are assigned fixed transition probabilities in terms of the inputs. A hierarchical graph is one where the lengths of the 
bonds and the transition probabilities obey scaling laws Q . 

The specific model we propose to study is based on a one-dimensional Lorentz lattice gas 0, H, the difference 
being that ours is a continuous time process where the separation between scatterers will be taken to satisfy a scaling 
law. The scattering probabilities depend on the direction of the particle, in analogy to the Lorentz lattice gas, with 
the further property that these probabilities change according to the index of the scatterer. Due to its connection to 
persistent random walks, we propose to refer to our model as a persistent hierarchical graph. In such a system, the 
evolution operator for phase space densities, the Frobenius-Perron operator, can be written explicitly and its spectral 
decomposition expressed in terms of the Pollicott-Ruelle resonances Sj . This in turn yields the expression for the 
survival probability : 



P{t)=Y,A i e« t , (3) 
j=o 

where the amplitudes Aj can be expressed in terms of the eigenstates associated to the corresponding resonances. The 
Pollicott-Ruelle resonance spectrum is located in the lower half-plane, Re Sj < 0. We note that the Perron- Frobcnius 
operator is here defined on a rigged Hilbert space, whose dimension is infinite pp[ . 

It is a general property that finite open graphs are fully chaotic systems. Indeed there is a gap empty of resonances, 
i. e. the closest resonance to the imaginary axis is real and isolated, so that it dominates the sum in Eq. (^|). This 
resonance is the escape rate. Thus in this case the survival probability decays exponentially as in Eq. (Q). 

On the contrary, in the semi-infinite open hierarchica l gra ph, because the lengths of the bonds and the transition 



probabilities scale in terms of some parameters (see Sec. [IB ), we will see that the decay of the survival probability is 
algebraic as in Eq. (0). Indeed the power law behavior can emerge in the limit of infinite graphs from the expression 
(|3|), because there is an accumulation of resonances going to zero and distributed with a particular density. In fact, 
if the amplitudes and the decay rates satisfy Aj — aa? and sj = — bft with a, (3, a and b some real functions of the 
parameters of the model (0 < a, j3 < 1), then, evaluating the sum in Eq. (||) by the steepest decent method we get 
the power law decay of Eq. (0) with 

Aj = aa3 \ tea 



We will show in Sec. IV that indeed these scaling behaviors for the spectrum Sj and for the amplitudes Aj hold for 
persistent hierarchical graphs. We further point out a connection between the parameters a and (3 and the scaling 
parameters of dynamical traps pp| , pT| : a is the spatial scaling parameter and 1/(3 the temporal one. We will 
comeback to this in the conclusions. 
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As already mentioned, every finite size approximation of a persistent hierarchical graph is a fully chaotic system. 
This provides means of making further comparisons between typical and fully chaotic systems. Dynamical quantities, 
such as the topological pressure (henceforth referred to as free energy) will be considered. 

Recent studies consider the question how the hierarchical structure and the dynamics of a typical Hamiltonian 
system shows up in quantum properties. It was shown by a semi-classical argument [22 that Eq. (pi) leads to fractal 



It was shown by a semi-classical argument 
conductance fluctuations on an energy scale larger than the mean level spacing. This fractal conductance fluctuations 
are decorated with peaks corresponding to isolated resonances at a small energy scale JT6|] , which are associated to 
"hierarchical" states and are distributed according to p(T) ~ 1/r, with T the width of the scattering resonance in the 
wavenumber plane, to be defined in Sec. 
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Quantum properties of graphs are not without their own interest. Spectral properties of closed quantum graphs 
have been considered in ]23[ | where quantum graphs were introduced for the first time as a model for quantum chaos. 
Other spectral properties were considered in 24 |2^, ^6|. Dynamics p7| , scattering [^8) and localization in infinite 
disordered graphs (29j have also been considered. On the other hand the classical dynamics has been studied in detail 
in PQ] . Here we will also study properties of the scattering resonances in a quantum realization of the persistent 
hierarchical graph. 

The plan of the paper is as follows. In Sec. [n| we discuss general properties of classical and quantu m g raphs and 



Secti 



introduce the persistent hierarchical graphs in some detail. The survival probability is defined in Sec 
give its expression in terms of the spectral decomposition of the evolution operator 
of the spectrum of the persistent hierarchical graph 



IV presents the calculation 
Some properties of the free energy are studied in Sec. [v| In Sec. 
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where we 



VI we turn to the quantum description of persistent hierarchical graphs and, in p articular, analyze the spectrum of 



scattering resonances. Finally conclusions and perspectives are drawn in Sec. VII 



II. HIERARCHICAL GRAPHS AND OUR MODEL 



A. General Survey 



A graph is a collection of B one-dimensional bonds, connected by vertices, where a particle moves freely. The 
position of the particle on the graph is described by a coordinate Xb- The index b refers to a particular bond and 
Xb to the position on that bond, with < Xb < h, where lb denotes the length of the corresponding bond. Here we 
consider oriented graphs where bonds have directions. Thus to each "physical" bond corresponds two oriented bonds, 
and therefore the number of oriented bonds is IB. 

In a quantum graph the dynamics of the particle on a bond is governed by the free Schrodinger equation. When 
the particle arrives to a vertex, a scattering process determines the probability amplitude abb' for being reflected or 
transmitted to the other connected bonds. 

These systems admit a classical limit which corresponds to a particle moving at constant velocity on the bonds and 
undergoing a conservative scattering process at the vertices (the classical limit of the quantum scattering process), 
determined by the probabilities Pbv = \&bb'\ 2 of being reflected or transmitted to other bonds. The classical dynamics 
is Markovian, i.e. there is no memory effect. 

Open or scattering graphs, have some infinite leads c, from where a particle can escape and never return. 

Hierarchical graphs are a particular class of graphs consisting of a self-similar collection of unit cells, which are 
topologically identical and whose characteristic lengths follow a given scaling law. In the classical case, the transition 
probabilities are taken to satisfy a scaling property, such as in a continuous Markov chain ||. A possible example 
is given by a random walk on a fractal support, such as the Sierpinsky gasket, where scattering probabilities change 
according to the level (with respect to the fractal structure) of the vertices. A simple such example is the persistent 
hierarchical graph we introduce below. Other examples are the quantum version of the chain model |l6| , or the Cayley 
tree for quantum conduction J3l], 



B. Persistent Hierarchical Graph 



The hierarchical graph that we consider is a semi-infinite one-dimensional lattice where the lengths of the bonds 
decay exponentially with the bond index. On this lattice, a random walker moves on the bonds with constant speed, so 
that the time between collisions becomes exponentially shorter as the walker moves deeper into the lattice. Moreover, 
at each vertex, the walker undergoes a random collision which reverses its direction with some probability q n , which 
depends on the index n of the scattering vertex, or keeps the direction of the walker unchanged with probability p n . 
We will label by n the (non-directed) bond between vertices n and n + 1. A directed bond b is either (n, +) or (n, — ). 
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FIG. 1: Schematic illustration of the possible transitions and their probabilities for a particle colliding with a scatterer. The 
parameter /i is here taken to be 1/2. There are a total of 5 bonds and 2 scattering leads in this example. 

In terms of those, the following transitions Pbv of going from b' to b are possible : 

y = * + \ y f b=(n + 1,+), with probability P (n+1)+)i („,+)= p n +i , ^ 

\ 6=(n,-), with probability P ( „ ,-),(«,+) = 9n+i ■ 

y = , n _\ > i b = {n-1,-), with probability P ( „_ 1 _ );(n _ } = p n , 

\ &=(n,+), with probability P(„ >+ ),(„ _) = q n ■ 

All the other probabilities are zero. 

The probabilities p n and q n are chosen so as to favor backscattering of particles as they move deeper into the lattice, 

Pn = Po£ n , , 7 j 
q n = 1 - p n , 

where < e < 1 is a fixed parameter. 

We study this infinite hierarchical graph by considering finite approximations made of N + 1 vertices with open 
boundary conditions at both ends. That is, infinite leads are connected to vertices 1 and N + 1 respectively from the 
left and right. For the length of the bonds l n we assume 

In = W l , l<n<N , (8) 

where Iq is an arbitrary length scale, which we will set to one, and /i, < /i < 1, is a (dimensionless) fixed parameter. 
We note that the total length of the lattice is Ln — Iq/j,(1 — — /j,), which is bounded by lofJ,/(l — n). A 

schematic representation of the system is presented in Fig. |l|. 

Since the particles move with constant speed, in the limit where N — > oo the particles will undergo exponentially 
more frequent collisions as they move to the right-most end of the lattice, and will be backscattered with exponentially 
increasing probability, hence, in practice, never reaching the right boundary. The escape is thus expected to be 
essentially due to exit from the left boundary for large enough TV. Figure || shows the results of a numerical simulation 
where the number of surviving particles is plotted vs. time. The power decay is apparent at long times, with an 
asymptotic exponent whose value agrees within a few percents with the value to be derived in Eq. (|3^). The 
configuration of the particles surviving after that time is displayed in Fig. |[ We note that the bonds are more or less 
equally populated at the exception of the left-most bonds, which are unpopulated. 



As discussed above Eq. (|l|) , the average (with respect to random initial conditions) of the ratio of the number of 
particles surviving after a time t to their initial number Nq defines the survival probability, which we would like to 
characterize in terms of the two parameters of the persistent hierarchical graph, namely e and /x. 
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FIG. 2: Fraction N(t)/No of particles remaining in the system after time t. The system has 25 bonds and iVo = 100,000 
particles are initially distributed at random positions (1 e. evenly with respect to the position on the line). The parameter 
values are e = 1/20 and /i = 9/10. Each particle is run for a maximal time 10 6 and escape times are recorded, which yields the 
fraction of particles of surviving particles vs. time. 

III. THE SURVIVAL PROBABILITY 

The statistical average of a physical observable A(xb) defined on the bonds of the graph is given by |3(J 

2B i fib 

< A >* -Er = (A\f*Po) , (9) 

6=1 lb J ° 

where pg denotes the initial probability density which evolves with the Frobenius-Pcrron operator P l to give, at time 
t, a density p(xb,t) at position Xb on the bond b. 

In particular, if we consider the observable A(xf ) ) = lb for the bonds that compose the finite part of the graph and 
A(x c ) — 0, for c a scattering lead, Eq. defines the survival probability, i. e. the probability of finding the particle 
in the interior of the system at a given time t, 

P(t) = J2 [" P(x b ,t)dx b . (10) 
b Jo 

We will henceforth reserve the notation A for this observable. One of our goals will be to show that this definition 
can indeed be decomposed as in Eq. (||). 

Since we are interested in the time evolution at long times, we may consider the spectral decomposition of P* to 
get an asymptotic expansion valid for t — > +oo of the form 

P(t) = (A\Ppo) = £m-> (%| Po ) + ■ • ■ (11) 

3 

as a sum of exponential functions 2 . Therefore Eq. (||) is obtained with 

A j = (A\* j )(* j \p ) . (12) 



2 Possible extra terms such as powers of the time multiplied by exponentials, t m exp(sjt), are not generic and may appear for particular 
values of the parameters of the system. See |B(| for a discussion of this point. 
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FIG. 3: Configuration of the fraction of particles surviving after t = 10 6 unit times. About 2/3 of the initial number of particles 
have escaped after the time considered. The particles distributed on the bounds with larger indices have essentially retained 
their initial positions. The escape occured from bound 1 only. 



The spectral decomposition used in Eq. ( |ll|) is fully determined by the solutions of the problem [|30| 

Q( s j)Xj = Xj , (13) 

which determine the Pollicott-Ruelle resonances Sj and the corresponding eigenstates Xj- Q is a 2B x 2B matrix with 
elements given by Qbb'(s) = Pbb<e~ slb ' ■ 

Explicit expressions for the scalar products in Eq. (|l^) were found in Q. For the right eigenstates, 



. «b Jo 



e Sj « A(x b ) dx b , 



and for the left eigenstates, 



E 6 - h»Xj[b"}*Xj[b 



(14) 



(15) 



Here \j [b] denotes the b component of the eigenstate Xj an d Xj denotes the complex conjugate of the left eigenvector 
ofQ(sj). 

If for the initial density we take po{x b ) = 1, for all bonds 6, and po{x c ) = for infinite leads c, that is a uniform 
distribution over the finite part of the graph, we have 



1 Y.b.b'XMxWKe 



.'\\ P s i l b' J- P -Sjh 



EblbXj(b)x*(b) 



(16) 



As we will show in Sec. |y|, the Pollicott-Ruelle resonances Sj are small so that we can expand the exponential terms 
in Eq. ([i"6|) and get, to first order, 



EbhxM Ebhx*[b] 
EbibxMxm 



(17) 
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IV. POLLICOTT-RUELLE RESONANCES 



According to Eq. (p_3[), the Pollicott-Rucllc resonances Sj are the roots of the following determinant, 

det [I - Q(s)] = . 

In order to write explicitly the matrix Q, we will order the states according to 



(1. 



1,- 2,+ 2,- 



N,+ N, — ) 



This way Q has the expression 
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Since this is a sparse matrix, it is rather straightforward to compute the determinant Eq. dig), 



N 



det [l-QOO] =Y[Si , 



(18) 



(19) 



(20) 



(21) 



where 



Si = 1 - qiqi+ie 



Si 



-2sli 



1 - qiq 2 e 



- 2 sh 



- 1 



U-l 



i > 2 , 



(22) 
(23) 



Owing to the product structure of Eq. (|T]), the zeros of Eq. (18) are the zeros of 5n- One can compute them 
numerically for any value of the parameters e and /i. It should be emphasized that the identification of zeros in Eqs. 
( pl| - |2"3| ) holds only for finite N . However, the numerical resolution of Eqs. (|22]-p3|) is limited to small N and, for the 
sake of proving Eq. (Q), a perturbative approach allows an analytic treatment. This is done in what follows. 



A. Perturbation Theory 



In order to set up a perturbation scheme, we choose e as our small parameter and note that the "unperturbed" 
system, e — 0, corresponds to the union of N non-interacting bonds. That is, the particles are oscillating back 
and forth on the same bond. The spectrum of this unperturbed system is the union of the s n , m — i^j^-, m € Z, 
n = 1, . . 
limit N 



, N. The resonances with s ^ are not degenerate and remain isolated under the perturbation, even in the 



oo. After the perturbation, each isolated resonance adds one contribution to Eq. (jllj), with an Sj whose 
real part is negative and 0(e). Therefore the states associated to them decay exponentially fast (with an oscillation 
on top). On the other hand, the resonance s = is the only degenerate unperturbed resonance, with a multiplicity 
N. As it will turn out, the perturbation acting on this resonance reduces the degeneracy by one unit at each order 
of the perturbation, the splitting being proportional to e J with j the order of the perturbation. Thus, in the limit 
N — * oo, the spectrum has an accumulation point at s = 0. Therefore these states cannot be considered isolated 
and their contribution to Eq. ( O ) must be accounted for separately from that of the isolated resonances because it 
becomes an integral in this limit. This integral accounts for the algebraic decay of the surviving probability in the 
long time limit. 

Let us discard the isolated resonances and consider only the resonances s nj Q. The eigenstates associated to the 
unperturbed system are solutions of the equation 



Q (0) (o)xi 0) =xl 0) 



(24) 
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where Q(°) is given by Eq. (EQ), in which e is set to zero. Explicit expressions for the eigenvectors are : 



(o) 
Xi = 



v/2 



1 







(o) 

X2 



1 





1 
1 





(0) 



1 

71 










1 

v i y 



(25) 



In order to implement the perturbation theory in powers of e, we first consider the right eigenvectors. The calculation 
transposes straightforwardly to the case of left eigenvectors. The perturbation theory closely resembles the standard 
perturbation theory for degenerate eigenvalues p4|. Let us consider linear combinations 



N 



(26) 



where the coefficients Ci are polynomials in e, and seek approximate solutions of the system 

Q(s)x = X , (27) 

where s is a polynomial in e and Q will be expanded to a given order in e. Writing Q = Q(°> + <5Q, we substitute 
Eq. @ into Eq. (§?]). Multiplying both sides of Eq. @ by X f \ i = h ■ ■ ■ , N, and using Eq. (E4), we obtain a 

system of N linear equations for the coefficients Cj, J2jLi Vi,j( s ) c j = 0, where V^j(s) = x- '' [Q(s) — Q^°H0)]Xj°' ) are 
the matrix elements of the perturbation operator 

/ -2+e- 3i i(g 1 +g 2 ) e-'2p 2 



vx(-) = i 



8,1 (91+92) 

Sil P2 





-2+e" 



^P2 
' 2 (92+93) 
' 2 P3 



-2+e" 



' 3 (93+94) 



V 





The values of s are found by solving the secular equation 

det[Vi(s)] 



\ 




-2+e- s! N( gjv+?JV + 1 ) / 



(28) 







(29) 



to the desired power in e. 

Expanding p„ and q n in powers of e, we can compute the corrections to the unperturbed solution. In fact, expanding 
Eq. (|29j) up to 0{e) we find that only one eigenstate, si,o, has negative real part, 



si.o 



' 2/i 



while up to this order, the others remain degenerate, 



Sn.o = 0{e 2 ) , n > 2 



The eigenvector corresponding to s^o is 



X(si) 



(o) 
Xi 



■0(e) 



(30) 



(31) 



(32) 



Hence the degeneracy remains to be lifted among the N — 1 remaining eigenmodes. We study now how the second 
order correction affects the degenerate state. We proceed in a similar manner as we did for the first order. The only 
difference is that now does not belong to the base of the degenerate subspace. Accordingly the perturbation 
operator in this subspace is represented by the matrix V2 = {Vi,j)2<i,j<Ni which is obtained from the matrix Vi by 
removing the first line and first column. Expanding the equation det[V2] = up to 0(e 2 ), we get a result similar 

to Eqs. (|30|-|32"|) with s 2 ,o = —poe 2 /2h + 0(e 3 ) and \i = X 2 °^ + C( e3 )- Proceeding, the effect of the perturbation at 
the third order in e must be studied among the remaining N — 2 degenerate states. By induction we thus have a 
whole hierarchy of roots, each corresponding to a different order in e and determined by secular equations involving 
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the corresponding perturbation operator that acts in the still degenerate subspace V„ = {Y%j)n<%j<N expanded up 
to 0(e n ). It is clear that, at any given order of the perturbation theory, the resonances which are not anymore part 
of the degenerate subspace will have further corrections to their values. However we do not need to take them into 
consideration since we are only interested in the leading contributions to every resonance of the spectrum. In fact we 
can prove the 

Proposition IV. 1 The N roots of Eq. ( |l8| ) can be approximated to order N in e by s^o, . . . , Sjv,0) where, for every 
1 < n < N, s ny o is the only root of order e" of the secular equation 

det[V„(s)] = , (33) 

with leading contribution 

Sn ,0 = -^f+O(e n+1 ). (34) 

The corresponding right-eigenvector Xi Eq. (p6|), has coefficients c n , . . . , cn which are the solutions of the linear 
system 

N 

^ Vj,k(s n ,o)ck = , j >n , (35) 

k—n 

Left-eigenvectors x have coefficients determined by 

N 

^ V jt k(s n ,o)cj = , k>n. (36) 

j=n 

We will not discuss the states associated to s n ,m (m 0)> which are exponentially decaying states and play a role 
only at the early stages of the dynamics. 

From Eqs. (|35|-^6|) it is easy to show that to the leading order we have 

Xn = Xn = Xi 0) + 0(s n ) . (37) 
B. Algebraic Decay 



According to Prop. [IV. 1| , the resonances s n ,o are 0{e n ) and therefore accumulate to s = as n becomes large, thus 
proving what we promised. Equation (|i~7| ) together with Eq. ( |37| ) and the expression of the unperturbed eigenvectors, 
Eq. (p5|), allow us to evaluate the leading contribution to Aj : Aj = 2/i- 7 +C(e), where we have substituted lj — Iq^ 
and Iq = 1. 

Turning back to Eq. (Q), we have shown that the decay is algebraic as in Eq. (|l|) with 

me/ mfJ,— 1 

We point out as a conclusion to this section that both parameters of the persistent hierarchical graph, e and [i, 
are necessary to grant the algebraic decay of the survival probability. This point will be further discussed in the 
conclusions. In what follows we will derive further properties of the persistent hierarchical graphs, first classical and 
then quantum. 

V. THERMODYNAMIC FORMALISM 

For the real time process we consider, the free energy (usually referred to as topological pressure) per unit time is 
defined in analogy to continuous time processes where the stretching factors are here replaced by the inverses of the 
transition probabilities at the vertices of the graph |30| : 

T{p)=lw. ^lnZ T ((3) , (39) 

1 — >oo 1 
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where the dynamical partition function 



Z T ((3)=Y / lPb b 1 ---Pi 



(40) 



is the sum over all trajectories of time length T (or equivalently length L = vT) of their respective probabilities raised 
to the power j3, with (3 > playing the role of an inverse temperature. We note that the length of the t raj ectories 
cannot be measured sharply because of the continuous time nature of the system. Rather the sum in Eq. @ should 
be understood as a sum over all trajectories whose lengths are within an interval T ± Ai, where At is fixed. In the 
infinite T limit, the value of At is irrelevant. 

Some properties of the free energy, are |llj : (i) T is a monotonically decreasing function of (3; (ii) T has a zero 
for some (3 = du, < dn < 1 (the strict inequality being due to the open boundaries), where (iii) du is the fractal 
dimension of the repeller with respect to a properly defined metric space (in the sense that the space of trajectories 
is a continuum where trapped trajectories form a subset with fractal dimension); (iv) for hyperbolic systems of one 
degree of freedom, — J 7 ' (da) is the value of the positive Lyapunov exponent on the space of trapped trajectories; (v) 
measures the rate of escape from the system; (vi) the difference between these last two quantities is the metric 
(Kolmogorov-Sinai) entropy on the repeller, Iiks = — ^'(dn)) and (vii) .F(O) = /itop is the topological entropy. 

As argued in pOfl , the free energy Eq. ( p9| ) can be obtained as the leading zero of the following zeta function 
C(s,/3) = det[l — Q^(s)], where Qp is identical to the matrix Q defined in Eq. (|20|), with the probabilities qi and Pi 
now raised to the power (3. Thus the free energy T((3) is the leading solution s of the expression Y[i Si — 0, where i 
takes values on the set of bonds and the Si are determined by the recurrence relation 



5i = 1 



P 



-2sL 



Si = l-Afe' 28 ' 1 



1+ - 



2t* 



Si-i 



- 1 



We 



01 - i - q 1 q 2 d 

prove the following asymptotic behaviors : 



(41) 
(42) 



Proposition V.l In the limit of large /?, the free energy T is linear in f3 with a coefficient exponentially small with 
respect to the number of bonds in the system, N , 



lim T{(3) = -/?i±£ (t 

fJ^oc 2 \/i 



N 



(43) 



The proof of this result follows by considering the probability of a particle bouncing off a given bond n for a time T. 
Let Wxin) be this probability. We have 

W T (n) = [ qn q n+1 } T ^ , 



1 - 



1+e f e 



o\ £ - 



Thus the ratio 



W T (n) 
W T (n - 1) 



-i T 



> 1 



(44) 



(45) 



is larger than 1, which implies that, as (3 — * oo, the free energy is dominated by particles bouncing off the last bond, 
i. e. F{(3) « lim T (l/T)ln[VK T (7V) /5 ]. Equation (||) follows. 

Proposition V.2 When f3 tends to zero, the free energy T has a limit independent of e, given by lim^o J~(f3) oc 
-^j- [1 + O(fi)] . Thus in the limit of large number of bonds N, the free energy has a singular limit, linijv^oo .F(O) = oo . 

This holds since the rate of creation of new trajectories per collision is the same at every site, while the rate of 
collision per unit time increases exponentially as the particles go deeper into the lattice. Hence the topological 
entropy is overwhelmingly dominated by particles bouncing off the right-most bond. In particular, the reason for its 
diverging with N — > oo is due to the existence of trajectories undergoing an infi nite number of collisions in a finite 
time, e. g. the trajectory ever moving to the right. Numerical evidence for Prop. V.2 is shown in Fig. ^[ 

We close this section with the observation that, in the infinite system limit, we expect the free energy to have a 
phase t rans ition at t he value j3 = 1. This should resul t fro m the asymptotic behaviors of the free energy discussed in 

oo, by Prop. |V.1| , the free energy is zero for every (3 > 1 since it is zero at 
it diverges as (3 — > 0, and thus must decrease steeply for < [3 < 1. Wc 
thus infer that matching the two curves at /3 = 1 results in a discontinuity of one of the derivatives of T . 



Props. |V_1| and [V^ In the limit N 
[3 — > oo. On the other hand, by Prop 
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FIG. 4: Topological entropy vs. N for three different values of /i : /i = 1 (solid line), 9/10 (dot-dashed line) and 8/10 (dashed 
line). The two dotted curves are proportional to l/fJ, N . 



VI. THE QUANTUM HIERARCHICAL GRAPH 



In this section we wish to explore the possibility of having a scaling relation for the widths of the resonances, as was 
suggested in Numerical analysis of the time evolution of quantum systems as considered in [^7], is possible 
for a finite syst em, but goes beyond our scope. Thus we consider a quantum system whose classical limit is the one 



defined in Sec. [IB. The quantum system is a linear chain with transition and reflection probability amplitudes 



cr (n+l,+),(n,+) — %/Pn+l 



0-(n ,-),(n,+) - V9«+l > ( 46 ) 
f(n-l,-),(n,-) = VP" ' 
cr (n,+),(n,-) = ^V^" ' 

with p n and (7„ defined as in Eq. Obb f = for all other possibilities. It is clear that at each vertex the scattering 
matrix is unitary and therefore the quantum problem is well defined. Moreover we have that Pw = \o~bb'\ 2 which 
shows |2^j that the classical limit of this quantum problem is indeed given by Eqs. (|^-||). 

The time evolution of a wave packet in an open system is controlled by the scattering resonances defined in 
the complex plane of wavenumbers k as the poles of the scattering matrix. For a quantum system, denoting the 
wavefunction by ip(xb,t), we can write the survival probability as : P QM (t) = J2b I dxb\ip(xb, t)\ 2 . Letting ip(xb,t) = 
X)r c r(zb)e~ lErt , where c r (xb) are determined by eigenstates of the evolution operator with complex eigenvalues E r , 
we can rewrite the survival probability as follows : 

P QM (i) = J2 CrCr , e *^-v)t e -(r P +r r ,)t/2 ? (47) 

r,r' 

where we have used the decomposition E r = kf, = e r — iT r /2, 

Since for short times the quantum evolution follows the classical one, it is interesting to study the distribution of 
scattering resonances in hierarchical graphs and look for manifestations of the algebraic decay in a quantum spectrum. 
Hufnagel et al. |l(| showed, in the framework of a quantum version of the chain model, that the distribution p of 
the width V = —4Reklmk of the quantum scattering resonances k satisfies p(T) ~ 1/F, V <C 1. Moreover, using an 
argument based on perturbation theory, they argued that the imaginary parts of the quantum scattering resonances 
satisfy a scaling relation. 

Given a scaling relation for the resonance widths, 1^ = /*, the widths distribution follows: 

P (T) = J2 spi - r) « J S(r r)dx ~ i (48) 
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This distribution has been associated to peaks that decorate fractal conductance fluctuations observed in energy 
scales larger than the mean level spacing. But since the width and height of the peaks in the conductance are 
determined by the imaginary part of the scattering resonances, the scaling behavior of the resonance widths is also 
contributing to the self similar, i. e. fractal, shape of the conductance. 

The scattering resonances are the zeros of the zeta function det[l — R(fc)] , with the matrix R(k) obtained by replacing 
s by ik and Pbv by aw in Eq. (^p|). As for the classical resonances, we can develop a perturbative approach in order 
to determine the quantum resonances and the corresponding eigenstates, R(k)<p — ef>. 

As opposed to the classical case, the zeroth order resonances are generally non-degenerate. Indeed a straightforward 
calculation of det[l — R(fc)] with e = yields the roots 



,(Q) _ (ZP+V* 



• P^- (49) 



Similarly to Eq. (E5j), the eigenvector corresponding to k niP is given by 



u n,p 



2n-l 2n 

(0 ... 1 (-l) P ... 0). (50) 



Unless fx = 1, the zeroth order quantum resonances are all isolated, whereas, for (j, = 1, kn. p = (2p + 1)tt/2 is 
independent of n, so that, for every different p, the resonances have an Nth order degeneracy. We consider the two 
different cases separately. 



A. Non-hierarchical Graph : fi = 1 



For (jl = 1, the situation is similar to Sec. IV. We can proceed by analogy and show the following 
Proposition VI. 1 For equally inter spaced scatterers, i. e. fj, = 1, the quantum resonances are given by 

fcn,p= (2p + 1) ^ -»^ + 0(e" +1 ), peZ. (/, = !) (51) 

Thus, for every integer p, the fc„,p satisfy a scaling law and we have an accumulation point at (2p + l)n/2. Since 
the inverse of the lifetime is T = — 4Rc k Im k and v = 2Re k is identified with the speed of the particle we have that 

F = V 23.p3 — a . 

A n,p 2 3 ' 

This result shows that asymptotically the classical and quantum lifetimes of the resonances with longest lifetime 
coincide in the graph with evenly inter spaced scatterers (/i = 1). This is in opposition to fully chaotic graphs where 
the strict inequality is satisfied 7 > T mini ]37j). However there is no algebraic decay in this case. 



B. Hierarchical Graph : /1 7^ 1 

The case \i ^ 1 is trickier. Indeed, one expects that the lowest order correction to k n p is 0(e n ), but it can only be 
determined in the perturbation theory provided we know the n — 1th order correction to the roots k n i y with n' < n, 
as well as the corresponding eigenvectors. In the remaining of this section, we will outline the derivation of the first 
order resonances and their eigenvectors and present in table | the results of a computation of k n ,p to fourth order of 
perturbation theory. 

In order to find the solution of R(k niP )<p niP = 4> n ,p, we will write again R(fc„ iP ) = R^(k^ p ) + SR(k n , p ). Upon 

expanding the eigenvectors 4> n ,p in terms of the basis spanned by the zeroth order eigenvectors (j)^ and (f)„\, 4> n , P — 

Sm=i D 9 =o.i c n,p,m,q^m, q , we will make use of the property that 4>m,q is an eigenvector of R (0) (fci°],), R (0) (fcl°p)(/>m^ = 

A- n ,p,m,q(f>m]q, with eigenvalue A„ iPim! g = (— l) 9 iexp[(— l) p+1 i7r/i m ~ n /2]. The first order correction to k„} p is the 

solution kn} p of the equation (pn^p 5R{kn]p + £k^}p)4>n p = 0, where 5R must be expanded to first order in e. The only 
non-zero first order corrections have n = 1, cf. Table ffl. The corresponding corrections to the zeroth order eigenvectors 
are given by 

1 (b^ T SR(k {0) +Fk {1) )ch i0) 
c nm =- ^ d ^P + n ' p) ^ P , K#M, (52) 



13 



K 

2 p. 

TV 

TV 

— 77 

7T_ 

— _ZI_ 
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7T 

2p 4 

— 7r 
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4 M 










0(e 2 ) 



ip [2-2e" r >'+Po + <;" r "po] 

8[l+ e *"M]p 
ip [-2+2e" r ^+p + i; i7r >'po] 



8[l + e"l'] / i 
[-l + e "/M ]po 
4[ 1 _|_ e ijrn] M 2 
[_l + e ^/M] 
' 4[l + e"'*'] M 2 



P()M 



" 12[l + e»' r M]2 (i 

. [l + ^ i ' r> '] 2 pg-3e"^pgM 
' 12[l + e" r M]2 fl 

[ l + e ^/M]2 [ _ 1 + e .x M]po + e i^/^ [1 + e i^ > , lp 2 

4[l + e "/M]2[ 1 + e ix M ] M 2 
[ l + e ^/M ] 2 [ _ 1 + e ^M]p + e i^/M [1 + e i^M ]p 2 

4[l + e"/f'] 2 [l + e i ' r ' J ]fi 2 



-1+e'' 



']PO 



' 4[l + e i7r /*']M 3 
[_l + e ^/M ]po 

4[l + e iT/M] M 3 





n 


P 


0(e*) 


1 
1 

2 
2 
3 
3 
4 
4 


1 


1 



1 



1 




,[l + e "'^]' 1 p*-2[-l+ e "^]p 2 [l+e 2 " r ^-2 e "^(-l+, 1 )] + e "'^p i ;[-2+ e "'M ( _2 +A1 )- M ], 1 


ld[l + e' J '^]'- s l i 

.[l + e"'^]3pjl+2[-l + e i ' r f]pi[l + e 2 " r ^-2 e i ' r "(-l + M)]- e " rA 'pg[2 + e""'(2+M)-M]M 

„.pg[-2n + (4-2n-poM)e i ' r/A '-(4-poM+3poM)e 2 " r/,1 +2Me 3i '" J ] 
lB[l+e*' r //*]3 M a 

,pg[-2n + (4-2n+3poM)e" r/ ' i -(4-pop-pop)>! 2l,r/ ' i +2Me 3 "''] 
16[l+e i,r /f»]3fi3 
[_l + e**>]j, 
4[l + e" r ' J ]n :i ) 

'4[l+e'"-C]/» 3 ) 
[_l + e>"/f']p 
4[l + e"/M] M 4 

[_l + e ™/M] po 

4[l + e i "/f'lM 4 



TABLE I: Quantum resonances up to fourth order in e. The zeroth order resonance fc„,j> is here written modulo 27r//i n . Notice 
that the second order correction to &2,p is purely real. The same holds for the third and fourth order corrections to kz iP . This 
suggests that the imaginary parts of k n ,p are 0(e 2n_1 ). 



which are different from zero for (m = 1, q = 1 — p) and (m = 2,q = 0, 1). One can proceed to higher orders along 
these lines. The results for k n>p are presented in Table [| up to fourth order. We note that our results suggest that the 
imaginary parts of the k ntP are C(e 2 " -1 ). Given that the distribution of the real parts of the k„, lP is rather uniform, 
this imply that the widths T scale identically to the imaginary parts of the scattering resonances. Hence Eq. 
seems to hold for this example. 



VII. CONCLUSIONS 



We have presented a simple model of an open hierarchical graph for which an analytical treatment of the algebraic 
decay of the survival probability is possible. The novelty of our approach lies on the successful application to the 
persistent hierarchical graph of a formalism originally developed in the framework of fully chaotic systems, where the 
survival probability decays exponentially. 

For the classical system, the computation of the survival probability was done using the spectral decomposition of 
the evolution operator. We showed that the algebraic decay relies in an essential way on the scaling properties of both 
the Pollicott-Ruelle resonances and their amplitudes. Using a pertubative approach we argued that the resonance 
spectrum has an accumulation point at the value zero, which is characterized by a scaling property in terms of powers 
of the expansion parameter. The structure of the corresponding eigenstates with respect to the length scales of the 
system yields the scaling of the amplitudes. 

This result must be contrasted to the observation of algebraic decay in the self-similar Markov chains |3], |l6| . 
Although the exponents are identical, the hierarchical graph is a dynamical process where randomness is involved only 
through the modelization of the collisions with scatterers, as opposed to self-similar Markov chains where transitions 
between states lack the spatial structure of our system. In the persistent hierarchical graph, the geometric role of the 
parameter [i is very clear, whereas in the self-similar Markov chains, fi represents an area which affects the transition 
probabilities between states. 
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As already pointed out in the introduction, the parameters e and \i define a hierarchical dynamical trap |2^, pl| , 
in the sense that /i is the ratio between successive length scales and /i/e the ratio between the corresponding staying 
times. Our result Eq. (38) is another instance of the relation of these parameters to the transport properties of the 
system, in this case the algebraic decay that characterizes the survival probability. It would be interesting to know 
what relation does this bear to the transport exponent of anomalous diffusion. 

Other aspects of the properties of the classical persistent hierarchical graph were studied through the application of 
the thermodynamic formalism. We computed the free energy (or topological pressure) per unit time in terms of the 
leading zero of a zeta function defined in analogy to discrete time systems. Different asymptotic regimes were studied. 
In particular, the topological entropy, which is the infinite temperature limit ((3 — > 0) of the free energy, increases 
exponentially with the number of bonds in the graph. In the limit of large number of bonds, the low temperature 
(,3 > 1) free energy tends to zero exponentially with respect to the ratio e/fi < 1. Moreover these results suggest 
that the free energy undergoes a phase transition at /3 = 1. 

For the quantum system, we used methods similar to the classical case and conjectured that the widths of the 
quantum scattering resonances follow a scaling law, in agreement with the numerically observed width distribution 
|p.6| . This argument was motivated by the computation of the resonances to the first few orders in perturbation 
theory. The limitation of this result, due to the complexity of the resolution of the quantum problem, illustrates the 
gap that separates the understandings of the classical and quantum approaches. The resolution of this question is 
open to future research by Bob Dorfman and others. 
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